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We use the recently developed typical medium dynamical cluster (TMDCA) approach [Ekuma et 
al, Phys. Rev. B 89, 081107 (2014)} to perform a detailed study of the Anderson localization 
transition in three dimensions for the Box, Gaussian, Lorentzian, and Binary disorder distributions, 
and benchmark them with exact numerical results. Utilizing the nonlocal hybridization function 
and the momentum resolved typical spectra to characterize the localization transition in three di¬ 
mensions, we demonstrate the importance of both spatial correlations and a typical environment for 
the proper characterization of the localization transition in all the disorder distributions studied. As 
a function of increasing cluster size, the TMDCA systematically recovers the re-entrance behavior 
of the mobility edge for disorder distributions with finite variance, obtaining the correct critical 
disorder strengths, and shows that the order parameter critical exponent for the Anderson localiza¬ 
tion transition is universal. The TMDCA is computationally efficient, requiring only a small cluster 
to obtain qualitative and quantitative data in good agreement with numerical exact results at a 
fraction of the computational cost. Our results demonstrate that the TMDCA provides a consistent 
and systematic description of the Anderson localization transition. 
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I. INTRODUCTION 

Disorder is ubiquitous in materials and can drastically 
alter their properties, in particular, their electronic struc¬ 
ture and transport properties, and even may induce elec¬ 
tron localization. This phenomenon is known as the An¬ 
derson metal-insulator-transitionjii 5 Here, the transition 
from a metal to an insulator is not characterized by the 
vanishing of the charge carrier density but by the can¬ 
cellation of the hybridization paths accompanying the 
quantum localization of the wave functions due to coher¬ 
ent backscattering from random impurities, deep-trapped 
states, etc. As a result, electrons that occupy such expo¬ 
nentially localized states are restricted to finite regions 
of space, and hence cannot contribute to transport. The 
Anderson insulator is gapless indicating that the single- 
particle excitations are essential in determining its physi¬ 
cal properties especially at low energies. While there has 
been significant progress in the quest to understand this 
phenomenon, a proper effective mean-field treatment is 
not yet fully developed. 

There have been numerous theoretical studies of dis¬ 
ordered electron systems employing computational tech¬ 
niques of varying complexity not limited to numerically 
exact methods including exact diagonalization, trans¬ 
fer matrix, and kernel polynomial methods*^ - — var¬ 
ious renormalization group techniques and mean- 
field theoriesJ^r— While numerical exact methods have 
been successfully used to study Anderson localization, 
they generally require the treatment of large clusters and 
the use of powerful computers. As a result, they are dif¬ 
ficult to extend to the treatment of interacting systems 
or chemically specific models. An alternative approach is 


offered by mean-field theories such as the coherent poten¬ 
tial approximation and its extensions : 23 ’ 24 ! 26 They map 
the lattice onto relatively small self-consistently embed¬ 
ded clusters. These methods have been successfully ex¬ 
tended to the treatment of interacting disordered systems 
and to chemically realistic models. Unfortunately, these 
methods have been woeful in the treatment of the Ander¬ 
son localization due mainly to the averaging procedure 
utilized and improvements in the environment describing 
the effective medium have been limited to single-sites. 

Due to the central role a mean-field theory that prop¬ 
erly describes the Anderson localization transition (ALT) 
may play for further progress in the study of electron lo¬ 
calization, there is a need to formulate such theory. Fur¬ 
thermore, a well-known long historical dichotomy exists 
between the mean-field results and the numerical data for 
the Anderson localization transition. Here, we demon¬ 
strate that the dichotomy can be reconciled by incorpo¬ 
rating spatial fluctuations in a typical environment into 
the mean-field theory thereby offering a solution to this 
long-standing problem and providing a systematic frame¬ 
work in the mean-field theory of the Anderson localiza¬ 
tion. 

The most commonly used self-consistent mean-field 
theory, in the study of disordered electron systems, 
is the coherent potential approximation (CPA)i 23 i 24 In 
the CPA, the original disordered lattice is replaced by 
an impurity placed in an averaged local (momentum- 
independent) effective medium. While the CPA suc¬ 
cessfully describes some one-particle properties, such as 
the density of states (DOS) in substitutional disordered 
alloys ) 23 ' 27 it fails to describe the Anderson localiza¬ 
tion transition. This failure stems, in part, from the 
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single-site nature of the CPA, as it is unable to cap¬ 
ture crucial multiple backscattering interference effects 
that lead to electron localization. Cluster extensions 
of the CPA, like the dynamical cluster approximation 
(DCA) 26 ’ 28 ’ 29 and molecular CPA2H, allow for the incor¬ 
poration of such nonlocal spatial correlations; however, 
they still fail to describe the Anderson localization tran¬ 
sition. The arithmetic averages of random one-particle 
quantities (e.g., density of states) calculated within such 
mean-field theories cannot distinguish between extended 
and localized states and are not critical at the Anderson 
transition^—— Hence, such average quantities cannot 
be used as an order parameter. This is the reason that 
most mean-field theories like the CPA-^i and its cluster 
extensions including the DCA, 28 - 29 ' 35 fail to provide a 
proper description of the Anderson localization in disor¬ 
dered systems. This failure is intrinsic to these theories 
as the algebraically averaged quantities, i.e., the aver¬ 
aged density of states, always favor the metallic state. 
This can be understood from the fact that in an infinite 
system of localized states, the average density of states 
is nothing but the global density of states, which is a 
smooth function of the disorder strength near the crit¬ 
ical point while the local density of states becomes dis¬ 
crete with a non-trivial system size dependence (see, e.g., 
Refsl25l [in. ,32i. Ifhil - IZol for a detailed discussion). Further, 
due to the large fluctuations in the local Green function, 
its typical value is far removed from the average one^i 
as such, the local average Green function, which char¬ 
acterizes these mean-field environments, does not have 
any information about the typical nature of the local or¬ 
der parameter needed to signal a phase transition. In 
sec. mm we will demonstrate the failures of finite clus¬ 
ter extensions of the CPA in characterizing the Anderson 
localization transition using the DCA. 

Finding a proper single-particle order parameter for 
the Anderson localization transition capable of distin¬ 
guishing between the localized and extended states is a 
major challenge in the study of disordered electronic sys¬ 
tems. In contrast to the arithmetic average, the geomet¬ 
rical averag e 25 ! 37 ’ 39 ! 42 ’ 43 gives a better approximation to 
the most probable value of the local density of states. Do- 
brosavljevic et al — developed the typical medium theory 
(TMT) to study disordered systems, where the typical 
density of states (TDOS), approximated using the geo¬ 
metrical averaging over disorder configurations, is used 
instead of the arithmetically averaged local density of 
states. They demonstrated that the TDOS vanishes con¬ 
tinuously as the strength of the disorder increases to¬ 
wards the critical point and it can be used as an effective 
mean-field order parameter for the Anderson localization 
transition. 

While the single-site TMT has been shown to be suc¬ 
cessful in describing localized electron states, it still has 
several drawbacks. In particular, it does not provide a 
proper description of the critical behavior of the Ander¬ 
son localization transition in three dimensions (3D). As 
a local approximation, the TMT neglects the effects of 


coherent backscattering and, as a result, the critical dis¬ 
order strength obtained is Wj MT = 1.65 instead of the 
numerically exact value W c ~ 2.1— i 10 i 44 ~ — for the box 
distribution (in our units). Also, the universal order pa¬ 
rameter critical exponent (which has also been called the 
typical order parameter exponent)/^ /3 obtained in the 
local TMT is /?™ T « 1.0 whereas its recently reported 
value is /3 ~ 1.67 j 46 ’ 50 Another crucial drawback of the 
single-site TMT in 3D is its inability to capture the re¬ 
entrance behavior of the mobility edge (energy separat¬ 
ing extended and localized electron states) demonstrated 
in exact numerical studie s 47 ’ 51 ’ 52 for the disorder distri¬ 
butions with finite variance: the box and Gaussian dis¬ 
order distributions. The TMT also underestimates the 
extended state regions in all the disorder distributions 
studied in this paper. 

The inadequacies of the single-site TMT can be reme¬ 
died by incorporating spatial (nonlocal) correlations by 
constructing its cluster extension. This can be achieved 
by using the DCA or molecular CPA schemes, which sys¬ 
tematically incorporate the missing nonlocal effects. 

In this paper, we show in detail a successful extension 
of the local TMT to its cluster version using ideas from 
the DCA. We demonstrate how the finite cluster exten¬ 
sion of the local TMT is able to systematically solve all 
the crucial drawbacks of the single-site TMT, indicating 
the necessity to include the missing nonlocal effects. One 
of the features of the Anderson localization transition in 
three dimensions is the non-sclf-averaging nature of the 
local quantities close to the localization transition, which 
obtain a highly skewed (log-normal) distribution. Hence 
special care must be taken in constructing a cluster ex¬ 
tension of the TMT. To avoid such self-averaging issues 
in the TDOS, we developed the typical medium dynam¬ 
ical cluster approximation (TMDCA) ji, which properly 
characterizes the Anderson localization transition in 3D 
and does not suffer from the self-averaging, by explic¬ 
itly separating the local part of the TDOS and treating 
it with a geometric average over disorder configurations. 
Hence, we are able to obtain a proper TDOS that charac¬ 
terizes the Anderson localization in 3D. We demonstrate 
the versatility of our method by applying it to the box, 
Gaussian, Lorentzian, and binary disorder distributions 
and benchmark it with numerically exact methods. 

The typical medium dynamical cluster approxima¬ 
tion scheme is demonstrated to be a systematic, self- 
consistent effective medium theory for characterizing 
electron localization. As a function of increasing clus¬ 
ter size, we demonstrate that the TMDCA achieves con¬ 
vergence of both the critical disorder strength and the 
trajectories of the mobility edge as a function of the clus¬ 
ter size. Furthermore, the TMDCA fulfills all the essen¬ 
tial requirements expected of a “successful” cluster the- 
ory, 28 ’ 53 We find that the TMDCA scheme is a systematic 
self-consistent effective medium theory to study Ander¬ 
son localization transition in three-dimensions, which i) 
recovers the original single-site TMT scheme at cluster 
size N c = 1; ii) recovers the DCA results at small dis- 
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order strength (when most states are metallic); iii) pro¬ 
vides a proper way to separate the energy scales such 
that the characteristic mobility edge behavior (for the 
disorder distributions with finite variance) is recovered; 
iv) captures the critical behavior of the Anderson local¬ 
ization transition with correct critical disorder strength 
W c and order parameter critical exponent /?, and pro¬ 
vides the correct description of the Anderson insulator 
at large disorder strength (when all states are localized); 
and v) fulfills all the essential requirements expected of 
a “successful” cluster theory<2 & £2 


The main problem addressed in this paper is how the 
mobility edge energies vary with disorder strength, their 
trajectories, and what happens to these trajectories in 
the proximity of the Anderson localization transition. 
Furthermore, since the DCA always becomes exact when 
N c —>• oo, the main role of the effective medium in ap¬ 
proaches based on the DCA, or its extensions, is to ac¬ 
celerate this convergence. For the ALT, we find that the 
effective medium formed from the average Green func¬ 
tion does not converge as N c becomes large, rather it is 
only able to describe the precursors to localization (cf. 
Section fill II) . However, we find that a number of effec¬ 
tive media based upon the geometrically averaged den¬ 
sity of states provide convergent results, i.e., an order 
parameter. So far, we find that the fastest convergence 
is provided by the TMDCA. 


The rest of this paper is organized as follows: following 
the introduction in this Section [H we present the model 
and describe the details of the formalism in SectionE] We 
present the results of our calculations in Section Hill The 
absence of localization in the DCA (shows only precur¬ 
sor to localization) is described in IHI 11 Detailed analysis 
of how we avoid self-averaging as the size of the cluster 
is increased is discussed in IIII 21 Then in IIII 31 we de¬ 
scribe how to treat the states close to criticality where 
the hybridization function vanishes leading to the devel¬ 
opment of poles (delta functions) in the imaginary part of 
the cluster-excluded Green function, &(K, w). A detailed 
analysis of our results for the box disorder distribution is 
presented in Sec. IIII Al Section HlI Bl presents our results 
for the binary disorder distribution. In Sec. IIII Cl we 
present our results for the Gaussian disorder distribution 
while Sec. IIII Dl shows the results of our computations for 
the Lorentzian distribution. In Sec. IIII El we discuss in 
detail the procedure for obtaining the critical parameters 
especially the order parameter critical exponent /? for the 
various disorder distributions, while in Sec. hue we ad¬ 
dress the discrepancy observed in the trajectories of the 
mobility edge energies at higher disorder strength. We 
summarize and discuss future directions in Sec. m In ap¬ 
pendices |E an d El we present a concise description of the 
developed transfer matrix and kernel polynomial meth¬ 
ods, respectively, used in benchmarking the TMDCA. 


II. METHOD 

We consider the Anderson model of non-interacting 
electrons subjected to a random potential. The Hamil¬ 
tonian (for spinless fermions) is given by 

H = - Y2 m 4 c . i + h - c •)+ ~ vhi- (i) 

<ij> i 

The first term is the kinetic energy operator due to hop- 
ping of electrons on a lattice. The operators cj(cj) cre¬ 
ate (annihilate) a quasiparticle on site i , and tij is the 
hopping matrix element between nearest neighbors (i,j). 
The second term is the disordered part parameterized 
by a local potential Vi, which is a random quantity dis¬ 
tributed according to some specified probability distri¬ 
bution P{Vi), rii = c\ci is the number operator, and /i 
is the chemical potential. We set 4t = 1 as the energy 
unit. In our analysis, we use different disorder distribu¬ 
tions. In particular, we consider the box (Bo), Gaussian 
(Ga), Lorentzian (Lo), and binary (Bi) distributions, re¬ 
spectively, with the corresponding distribution functions 
P{Vi) 


P B o{Vi) 

P Ga (Vi ) 

P Lo {Vi) 

P Bi (Vi) 


2W b , 


-Q(W B o-\Vi\), 


2 *W% a 

W Lo 


0 -3V 2 /(2W 2 a ) 


<v 2 + w 2 Lo y 

c a S(V i -W A ) + c b 6(V i -W B ), 


(2a) 

(2b) 

(2c) 

(2d) 


where 0(a:) is the step function, c a is the concentration 
of the host A atom, Cf, = 1 — c Q is the concentration 
of the impurity B atom, and the strength of the dis¬ 
order in units of 4 1 is parameterized by W (W a and 
W B for the binary alloy model). We have scaled W Ga 
such that the second moments of the P Bo and Pgo agree 
with each other (i.e., we set the variance of the Gaus¬ 
sian distribution equal to that of the box distribution: 
a 2 = W 2 / 3) in the event that Wb 0 =Wg 0 to enable 
comparison. Since the Lorentzian distribution lacks a 
second moment, the disorder values cannot be directly 
compared with that of either the box or Gaussian distri¬ 
butions. We introduce a shorthand notation for disorder 
averaging: (...) = /dV)P(V))(...). 

To solve the Hamiltonian ©, different methods will 
be used including the DCA^ the cluster typical medium 
theory^ and the recently developed typical medium dy¬ 
namical cluster approximation (TMDCA) X We will com¬ 
pare these results to those obtained from numerical meth¬ 
ods like the kernel polynomial method (KPM)^^ - — and 
the transfer matrix method (TMM)4*i2c— 

The TMDCA utilizes the self-consistent framework of 
the standard dynamical cluster approximation^ with the 
important usage of an environment defined by a typi¬ 
cal nonlocal hybridization function. In particular, the 
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TMDCA maps the given disordered lattice system onto 
a finite cluster which is embedded in an effective self- 
consistent typical medium. Note that unlike the usual 
DCA scheme, where the effective medium is constructed 
via algebraic averaging over disorder configurations, the 
TMDCA scheme uses geometric averaging. By mapping 
a d-dimensional lattice to a finite small cluster contain¬ 
ing N c = Ld sites, where L c is the linear dimension 
of the cluster, we dramatically reduce the computation 
effort^ Unlike the single-site methods commonly used 
to study disordered systems, such as the coherent po¬ 
tential approximation (CPA) 23 i 34 or the local TMT^ 
the TMDCA ensures that nonlocal spatial fluctuations, 
neglected in single-site approaches, are systematically in¬ 
corporated as the cluster size N c increases. Short length 
scale correlations are treated exactly inside the cluster, 
while the long length scale correlations are treated within 
the typical medium. 


Algorithm: The details of the TMDCA formalism are 
described below. The nonlocal (K-dependent) disorder 


average cluster density of states is given as 
p c avg (K,u>) = (p c (K,uj,V)) = -^(QG C (K,K, W ,U)), 

(3) 

where the superscript ‘c’ denotes cluster and (...) is the 
disorder average. For a single-site N c = 1, we recover 
the CPA. The K-dependent cluster Green function is ob¬ 
tained from the site dependent Green function G c (i,j,uj ) 
via the Fourier transform: 

G C (K,K, W ) = i-^ e ‘ K ( R *- R ,)G c (*,j». (4) 

Within the TMDCA, for each cluster configuration, 
we first obtain p c (K,w) = — 3G C (K, K, u>)/n. It can 
be shown via the Lehmann representation ; 54 : 55 that 
P c { K, w) > 0 for each K, u>, and disorder configuration. 

As mentioned above, in the TMDCA, the local part of 
the cluster-momentum-resolved typical density of states 
is separated and treated with geometrical averaging over 
the disorder configurations, to avoid self-averaging as the 
cluster size increases. The obtained cluster typical spec¬ 
tra are given by 


local TDOS 


Nc 


Pt 


yp\ 


K,w) = exp ( (Inp-(w, V)) 


N, 


p c (K,to,V) 

b £i Pi^’ v ) 


nonlocal 


( 5 ) 


From Eq.[5l the disorder averaged typical cluster Green 
function is obtained using the Hilbert transform 


ori, r o (K , w) set to a small imaginary number may 
serve as the starting point. 



dui 


,P c t VP { K V) 


( 6 ) 


A schematic TMDCA self-consistency is shown in Fig. |T| 
The TMDCA iterative procedure is outlined below. 


2. After setting up the cluster problem, we calculate 
the cluster-excluded Green function Q( K,w) as 


1. The TMDCA iterative procedure begins by propos¬ 
ing an initial guess for the hybridization function 
r o (K , w), where the subscript ‘o’ denotes old. The 
choice of the starting guess for the hybridization 
function may be based on a priori knowledge, i.e., 
having information about the self-energy E(K,w) 
and cluster Green function G c (K,w), r o (K,w) can 
be calculated as 

F 0 (K, w) = w - e(K) +p- E(K, ui) - 1/G C (K, tv), (7) 
where e(K) = N c /N E e ( K + k) is the coarse- 

k 

grained bare dispersion with k summed over N/N c 
momenta inside the cell centered at the cluster mo¬ 
mentum However, if nothing is known a pri- 


Solve Cluster 
Problem 


G C (K) 


FT:g(K)->g(n,m) 

£(iC)=(io-r„-e t +gb 

Hi_£*_ •__ 3e_ : 


r„(K) = r a (K) + t[(G c (K)) -(G(K)n 
typ 


+r o (K)-e(K + k)+e k ] 


E P (K)=exp( i:^ ln(p ' C)>x( i^Xp r > 



p c (K,U 3 )* 
G c (K)= f do)'-®- 

typ 


FIG. 1. (Color online) The self-consistent loop of the 
TMDCA. 
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S(K, W ) = (a;-r o (K,a;)-e(K) + M )~ 1 . ( 8 ) 

Since the cluster problem is solved in real 
space, we then Fourier transform Cy(K,w): 

Qn,m — 6Xp(iK • (t n r* m )). 

3. Next, to solve the cluster problem, we stochasti¬ 
cally generate random configurations of the disor¬ 
der potential V, and calculate the corresponding 
cluster Green function as 

G C {V) = (g- 1 - v r 1 . (9) 

This is Fourier transformed to G c (K,K,w) 
to obtain the cluster density of states 
p c { K,w) = -i»G c (K,K,w). ' The typical 

cluster spectra is then calculated via geometric 
averaging using Eq. 0 Then, we calculate the 
disorder averaged, typical cluster Green function 
G c typ ( K , u>) via the Hilbert transform using Eq. 0 
We note the advantage of the stochastic sampling 
of the disorder configurations. Here, each of the 
disorder configurations is statistically independent 
of the others. Thus, for example, for the binary 
disorder distribution, instead of enumerating all 
configurations, which scales as 2 Nc , we do a 
stochastic sampling of the disorder configurations. 
This greatly reduces the computational cost, at 
the expense of a small sampling error, enabling us 
to study larger clusters. We also enforce all of the 
cluster translational and point group symmetries, 
effectively generating more configurations. With 
this, the number of disorder realizations needed to 
obtain a converged solution falls with increasing 
cluster size. For a typical 64 site cluster, with 
box disorder, about 500 disorder realizations are 
enough to produce high-quality data. The code 
scales like A(N C )N j) due to the matrix inversion in 
Eq.0 but the prefactor A(N C ), also depends on N c 
since fewer self-consistency iterations and disorder 
configurations are needed for larger clusters. 
Hence, A(N C ) falls with increasing cluster size. 

4. After solving the cluster problem, we use the typi¬ 
cal cluster Green function, G% yp (K, w), to calculate 
the coarse-grained cluster Green function of the lat¬ 
tice G(K,w) as 

G(K w ) = /'_ N S (K,e)de _ 

1 J J (Gt yp (K, w )) _1 + F(K,w) — e + e(K) ’ 

( 10 ) 

where N§ (K, e) is the bare partial density of states. 

5. We then close our self-consistency loop by updating 
the new hybridization function using linear mixing 

Fn(K, W ) = r o (K, W ) + i[{G c typ { K,ce ))- 1 - (G(K, a;))" 1 ] 

( 11 ) 


where the subscripts ‘n’ and ‘o’ denote new and 
old, respectively. The mixing parameter £ > 0 con¬ 
trols the ratio of the new and old T(K, ui) entering 
the next iteration. For very small £, convergence 
may be slowed down unnecessarily, while for very 
large £, oscillations about the self-consistent solu¬ 
tion may occur. Instead of linear mixing, the con¬ 
vergence of the computations can be improved by 
using the Broyden method^ 

6 . We repeat the above procedure until the hy¬ 
bridization function converges to the desired ac¬ 
curacy, r o (K,w) = r„(K,w). When this hap¬ 
pens, the Green functions are also converged, 
G(K,w) = Gt yp (K.,uj) within the computational 
error. 

We note that our formalism preserves causality just as 
the DCA,— since all the Green functions are causal, both 
the average density of states (ADOS) and the TDOS cal¬ 
culated from them are positive definite. Also, we observe 
that as N c increases, our method systematically interpo¬ 
lates between the local TMT and the exact result. 

For reproducibility, we specify in Table U the clus¬ 
ter geometries and other important parameters of the 
clusters used in our computations. The parame¬ 
ters of Table HI include the lattice vectors ( 01 , 02 , 03 ), 
and the cubicity (C)j^ The cubicity is given as 
G = max(ci,cj" 1 ) x max(c 2 , dj" 1 ), where Ci = 3 1//2 //d 
and C 2 = 2 1 ^ 2 l/f are cluster parameters defined by the 
geometric mean of the lengths of the four body diagonals 
of the cluster, d = (did^dadi) 1 ' 4 ', the six-face diagonals, 
/ = {hhhUhh) 1/fi , and the edges, l = {lihh ) 1 / 3 
C = 1 is for a perfect cube, and G > 1 otherwise. Fol¬ 
lowing this criterion, clusters N c = 1, 64, 125, and 216 
are perfect cubes. 

TABLE I. Three-dimensional (3D) cluster geometries utilized 
in our calculations. The at denote the cluster lattice vectors 
and C is the cubicity. 


N c 

a,i 




03 


C 

1 

(i, 0 , 

0 ) 

( 0 , 1 , 

0 ) 

(0, 0 , 

1 ) 

1.000 

38 

(1, 2, 

3) 

(3, -1, 

-2) 

(2, -2, 

2 ) 

1.087 

64 

(4, 0, 

0 ) 

(0, 4, 

0 ) 

(0, 0 , 

4) 

1.000 

125 

(5, 0, 

0 ) 

(0, 5, 

0 ) 

(0, 0 , 

5) 

1.000 

216 

(6, 0, 

0 ) 

(0, 6, 

0 ) 

(0, 0 , 

6 ) 

1.000 


III. RESULTS AND DISCUSSION 

Before presenting our main results in detail, we will 
first review the characteristics of the DCA. Despite its ad¬ 
vantage over the CPA, it shares the same behavior with 
the CPA in that it is unable to detect the localization 
transition in a disordered electron system^ We will also 
elaborate on the details of how self-averaging is avoided 








6 


in the TMDCA as the cluster size is increased. This 
becomes imperative since in the cluster, self-averaging 
will ultimately destroy our ability to detect the Ander¬ 
son localization transition. Both the inability of the DCA 
to capture the Anderson localization transition and how 
self-averaging is avoided in the TMDCA will be demon¬ 
strated using the box disorder distribution. 

1. Absence of Localization in the DCA 

The dynamical cluster approximation (DCA), unlike 
the coherent potential approximation (CPA) incorpo¬ 
rates nonlocal spatial correlations systematically as the 
size of the cluster is increased. While spatial corre¬ 
lations are an important ingredient in the localization 
transition in disordered electron systems, the DCA ef¬ 
fective medium is characterized by arithmetic averaging 
over the disorder configurations. As explained above 
and will be demonstrated below, even a typically de¬ 
fined medium without a proper treatment of the typi¬ 
cal density of states (the local part of the typical den¬ 
sity of states needs to be separated and treated explic¬ 
itly using geometrical averaging), reduces to the DCA 



FIG. 2. (Color online). Left panels: The average den¬ 
sity of states calculated using the DCA for small and 
large disorder strengths^. Right panels: The proba¬ 
bility of an electron remaining on a site after time t, 
P(t) = {\G(l, l, t)| 2 ). The insets in the plots on the 
right are the probability of an electron remaining on a 
site for all time P(t —> oo ,rf) = lim (\G(l,l,t)\ 2 ) = 

£—>•00 ' ' 

= lirn tL de (\G(l, l, t + ir/)\ 2 ') for the N c = 1 (black) and 

64 (red), respectively. The ADOS is not critical at the An¬ 
derson transition. The non-critical behavior in the DCA can 
also be inferred from P(t), which decays rather fast instead of 
remaining constant, and in that P(t,r]) extrapolates to zero 
instead of a finite value when many of the states are localized. 
The DCA only shows a precursor to the Anderson localization 
as manifested in P(t) for large disorder strength for the finite 
cluster sizes. 


for large clusters. In Fig. [ 2 ] (left panel), we show the 
ADOS at small and large disorder strengths for the 
single-site (N c = 1 ) and finite cluster (N c = 64 ) cal¬ 
culations. We also show in Fig. [ 2 ] (right panel) the prob¬ 
ability P{t) = (|G(£, l, f)| 2 ) of an electron remaining at 
a site l at long time t calculated within the single-site 
and finite cluster DCA— scheme while the insets de¬ 
pict the probability of an electron remaining on a site 
l for all time: P(t —I 00,77) = lim (\G(l,l,t)\ 2 ) = 

t —>OO N ' 

= £/“ 0 de(|G(Z,Z,e + w?)| 2 ). 

P(00,77) is expected — to be nonzero for any fraction 
of the localized states in the spectrum of the eigenstates 
of the disordered system . 28 ! 53 : 58 : 59 As is evident from the 
inset of the right panel of Fig. [ 2 J the plot of P(oo,rj) 
versus 77, P(00,77) extrapolates to zero even very close to 
the critical disorder strength for both N c = 1 and 64 , 
respectively. 

Also observe from Fig. [2] (left panel), the ADOS is not 
critical at the Anderson transition. This is manifested 
in the P(t) plot since for a localized state, P(t) is ex¬ 
pected to be finite. For N c = 1 , P(t) falls quickly with 
time regardless of the disorder strength. However, for 
N c = 64 , the electrons remain localized for longer times 
as the disorder strength is increased. This can be under¬ 
stood by noting that each site on the cluster is coupled 
to a non-interacting translationally invariant host into 
which electrons can escape. Hence, if a finite density of 
states exists at some energy, the corresponding states are 
guaranteed to be extended unless the hybridization rate 
between the host and the cluster vanishes. Indeed, this 
is the case in the DCA. The imaginary part of the inte¬ 
grated hybridization (escape) rate (—3 f T(K, uj)dK.duj) 
between the cluster and the host as a function of disorder 
strength (W) remains constant regardless of the strength 
of the disorder (cf. insets of Figs. [6] [H] and [T6l for dif¬ 
ferent disorder distributions). It does not go to zero as 
needed for a localization transition in the typical medium 
context. Thus, the DCA is only able to capture the pre¬ 
cursor to the Anderson localization as the cluster size, 
but not the disorder strength increases. 


2. Avoiding Self-averaging 

The averaging procedure used to calculate the typi¬ 
cal spectra is not unique. As noted above, our initial 
attempt to formulate a cluster version of the TMT re¬ 
produced the expected behavior in one and two dimen¬ 
sions as N c is increased. 33 However, in three dimensions, 
applying the algorithm directly will lead to an effective 
self-averaging for large clusters. This is due to the fact 
that close to criticality, there exist distinct localized and 
extended states above and below the localization edge 
given by the TDOS with an energy scale difference that 
can span an order of magnitude. These energy scales 
need to be treated differently. This can be seen by in¬ 
vestigating the spectra where the local part of the TDOS 
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is not explicitly separated and treated with a geometric 
averaging over disorder realizations: 


Pt yp ( K ^) = exp (Inp c (K, u, Vi)). 


In forming the Fourier transform 


( 12 ) 


p c {K,u,Vi) 



7 r 



53 e iK < x - x ' ) G c (X,X',uj,Vi) 

X,X' 

(13) 


we average over the cluster coordinates X and X' , in¬ 
cluding the local part, X = X' . So, the local DOS 
is first averaged over the cluster sites and then Fourier 
transformed making the local part of p c ( K,w). Hence, 
for large clusters, this reduces to linear averaging of 
the local part instead of geometrical averaging. As a 
consequence, the host Green function constructed from 
Pt yp ( K,w) of Eq. [12] is unaware of the TDOS and thus, 
it is unable to distinguish between the energies above 
and below the localization edge. To avoid such self¬ 
averaging in the TDOS, we proposed the Typical Medium 
DCA (TMDCA) method^ Here, the cluster-momentum- 
resolved typical density of states (TDOS) for each K 
is split into local and nonlocal parts. The local part 
is treated with geometrical averaging over disorder con¬ 
figurations, while the nonlocal part is treated with an 
algebraic or geometric averaging over the disorder con¬ 
figuration. 

To do this, we have utilized two schemes. The first 
scheme is what we call linear-log procedure, which is used 
in this study. Here, we treat the local part with a geomet¬ 
rical averaging while the nonlocal part is approximated 
algebraically using linear averaging as 


Ptypi. K ,u) = exp 53 (In Pi(u, x 

iT"; Vi U )- (l4) 

W c ^ fob V i) / 

The second scheme is what we call the log-log proce¬ 
dure, which again involves the treatment of the local part 
with geometrical averaging; however, the nonlocal part is 
treated with a log averaging as 


Pt yP (, K >u) = ex P 



* ” c 

53 (InptfaVl)) 


x 


exp 



P c { K,ce,H z ) \\ 
W c £i Pi ( w > V i) 1) 


(15) 


It is imperative to note that while there are dif¬ 
ferent behaviors of the two methods around the re¬ 
entrance region (cf. Fig. [3]), both procedures systemati¬ 
cally converge to the same critical disorder strength e.g., 
W^ c - 12 ~ 2.1 ± 0.01 for the box disorder distribution. 


However, the former ( linear-log procedure) is generally 
more robust than the log-log method. The latter is char¬ 
acterized by slower convergence around the re-entrance 
region, requiring far larger cluster sizes before the conver¬ 
gence of the re-entrance region is achieved in comparison 
to the e.g., the transfer matrix method (TMM) results. 
It may also not be adequate to study localization phe¬ 
nomena in realistic material applications, since it is not 
obviously clear how a geometrical averaging of the off- 
diagonal components of the spectral density, which are 
not positive definite, will be done. The comparison of 
the phase diagram obtained using the two procedures: 
log-log and linear-log formalisms is shown in Fig. [3] As 
is evident from the figure, the two new schemes converge 
to the same critical disorder strength but behave differ¬ 
ently around the re-entrance region. While our origi¬ 
nal formulation (cluster typical medium theory (CTMT)) 
will eventually converge to a disorder strength far greater 
than W c , a further remark is that the re-entrance trajec¬ 
tory of the mobility edge is totally missed as a conse¬ 
quence of self-averaging in the cluster. 



FIG. 3. (Color online). A comparison of the phase di¬ 
agrams of the Anderson localization transition in 3D ob¬ 
tained from different cluster approximations with N c = 38 
using the CTMT and the TMDCA (linear-log and log-log) 
schemes. Observe that in the CTMT, as a consequence of 
self-averaging, the higher disorder behaviors that are captured 
in our TMDCA are totally missed and the critical disorder 
strength is also severely over-estimated. 


We note that in both the linear-log and the log-log 
procedures, at small N c , about 100 self-consistent itera¬ 
tions are required to achieve convergence, while, for rel¬ 
atively large N c , far fewer iterations are required. The 
convergence criterion in both limits is achieved when the 
TDOS [ui = 0) does not fluctuate anymore with iteration 
number within the error bars. 

Finally, we note that many other definitions of the typ¬ 
ical medium which avoid self-averaging are possible, in- 








eluding the use of only the local part of Eq. [Ml i.e., 


Pt yp ( K ,u) = exp ^ (In Pi(w, Vi)) 


(16) 


However, this method was rejected since it does not meet 
most of the criteria discussed in the introduction. In this 
case, this formalism does not recover the DCA in the 
weak disorder limit. 


3. The Pole Procedure 


Close to criticality, the hybridization rate of the states 
at the top and bottom ((0,0, 0) and ( 7 r, n, 7 r)) of the bands 
tends to zero leading to the development of poles in the 
cluster excluded Green function. Here, we present in 
detail how to deal with such poles that emerge on the 
real frequency axis as the critical disorder strength is 
approached. 

When 3r i j /p (K,cc) becomes very small (i.e., 
QEi yp (K, w) —> 0), the imaginary part of the cluster- 
excluded Green function, £(K,w), becomes a series of 
delta functions. To see this, we note that 

G(K, w) = (w - T typ (K, w ) - e(K) + m)" 1 (17) 

= V(u> — a/) -1 — iTr5(u} — cj'), 

where u/ = e(K) — p + 5?T typ (K,u;) and "V‘‘ denotes 
the principle value. Evidently from Eq. 1171 the poles 
cannot be represented in the conventional way as a list 
of frequencies on the computer with a finite frequency 
resolution dui. Such difficulty can be avoided by replacing 
G(K,lu) for each of the K—cells where 3r typ (K,u;) is 
vanishing with 


G(K,u>) 



uj = oj 1 

OJ ^ U)’. 


(18) 


We refer to this formulation as the explicit “pole- 
procedure". With this procedure, the singularity in 
Q(K,oj) can be properly captured. An added difficulty is 
that for a given N c , as the W c is approached, 3r ipp (K, oj) 
for individual cells goes to zero at different rates. Hence, 
we have to determine which of these cells need to be 
treated with the explicit “pole-procedure”. We choose 
the criterion that for any cell, if (— 1/7r) x 3r tpp (K, oj') 
< a x doj', then, we apply the pole procedure to such cells. 
Here, a > 1 is a parameter that measures the minimum 
number of pixels needed to represent a pole approaching 
the real frequency axis. Our numerical experience shows 
that such a criterion works nicely while spurious results 
are obtained otherwise. 


A. Box Disorder Distribution 

To demonstrate that the typical and not the average 
DOS can serve as a proper order parameter for the An¬ 
derson localization transition, we start the discussion of 


our results by comparing the algebraically averaged DOS 
(ADOS) calculated using DCA and the TDOS obtained 
from a single-site TMT (N c = 1) and finite clusters 
TMDCA (N c = 64 and 125) at various disorder strengths 
W for the box disorder distribution lEa. l2al) . As shown 
in Fig. [4] the ADOS remains finite while the TDOS for 
both TMT and TMDCA continually gets suppressed as 
the critical disorder strength is approached. Moreover, 
one observes a crucial difference between the single-site 
TMT (N c = 1) and TMDCA finite clusters of N c = 64 
and 125. In the former, the mobility edge (for extended 
states TDOS is finite) defined by the boundary of the 
TDOS (indicated by arrows) always gets narrower with 
increasing disorder strength W, while in the latter, as 
a function of disorder strength, the mobility edge first 
expands and then decreases, hence giving rise to the re¬ 
entrance behavior, which is completely missing in the 
single-site TMT. Observe also the quick convergence with 
the clusters size at ui = 0 for the finite clusters N c = 64 
and 125. The implications of this will be discussed fur- 


N =1 N =64 N =125 

c c c 



FIG. 4. (Color online). The evolution of the ADOS and 
TDOS at various disorder strengths W for the single-site 
TMT and the TMDCA with cluster size N c = 64 and 125. At 
low disorder, where all the states are metallic, the shape of 
TDOS is the same as that of the ADOS. As W increases, in 
the case of single-site TMT, the TDOS gets suppressed and 
the mobility edge moves towards w = 0 monotonically. In the 
TMDCA, the TDOS is also suppressed, but the mobility edge 
first moves to higher energy,and only with a further increase 
of W > 1.8, it starts moving towards the band center, indicat¬ 
ing that the TMDCA can successfully capture the re-entrance 
behavior. Arrows indicate the position of the mobility edge, 
which separates the extended electronic states from the local¬ 
ized ones and the colored region indicates the TDOS. 
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FIG. 5. (Color online). Comparison of the average (typical) 
density of states calculated with the DCA (TMDCA) and the 
kernel polynomial method (KPM) for the box disorder distri¬ 
bution at various strengths for the cluster size iV c =64. The 
kernel polynomial method used 4096 moments on a 48 3 cubic 
lattice, and 1000 independent realizations generated with 32 
sites randomly sampled from each realization. 


ther below and explored with respect to what happens 
to the trajectories of the mobility edge as the size of the 
cluster increases. 

To benchmark our TMDCA formalism with another 
numerical technique, we utilize the kernel polynomial 
methodi^ We show in Fig.[5]a plot of the TDOS (ADOS) 
(N c = 64) calculated using the TMDCA (DCA) as com¬ 
pared to the TDOS (ADOS) obtained using the kernel 
polynomial method (KPM). As it is evident from the 
plots, even though there is a qualitative agreement be¬ 
tween the two methods, there are subtle deviations es¬ 
pecially in the TDOS. This deviation can be attributed 
to finite lattice effects and the effective broadening due 
to the finite order expansion used in the KPM. Over¬ 
all, the agreement is a manifestation of the ability of our 
TMDCA formalism^ to accurately characterize the An¬ 
derson localization transition in systems with a uniform 
disorder distribution even with relatively small system 
sizes, as compared to the large lattice systems that need 
to be simulated in the kernel polynomial method for ac¬ 
curate results to be obtained. 

Next, we consider the evolution of the critical disorder 
strength W c with the cluster size. Figure [6] shows the 
local TDOS(w = 0) at the band center as a function of 
disorder strength W for several cluster sizes: N c = 1, 64, 
and 125. The critical disorder strength W c is defined 
by the vanishing of the TDOS(w = 0). The inset is the 
imaginary part of the integrated hybridization function 
which shares the same property as the TDOS since both 
vanish at the same disorder strength, while the DCA re¬ 
sult remains finite, independent of the disorder strength, 



FIG. 6. (Color online). The TDOS(u; = 0) vs. dis¬ 
order strength W at different cluster sizes N c — 1,64, 
and 125 for the uniform (box) disorder distribution. The 
TDOS(.R = 0,a; = 0) vanishes at the critical disor¬ 
der strength W c - At N c = 1, lFj Vc=1 ss 1.65. As 
the cluster size increases, W c systematically increases with 
W / J Vc » 12 ~ 2.10 ± 0.10, showing a quick convergence with 
the cluster size. The inset shows the integrated hybridiza¬ 
tion function (—Ss f T(K, co)dKduj) as a function of disorder 
strength W. Observe that — $y f T(K, uj)dKduJ vanishes at 
the same disorder strength as the TDOS. The dashed lines 
are the —$y f r(K,w)dKdcj from the DCA. Observe that it 
is a constant regardless of the disorder strength and clus¬ 
ter size. This shows that the DCA, even though it incor¬ 
porates spatial correlations, does not describe the Anderson 
localization transition. Moreover, near the critical region the 
TDOS(.R = 0,o i = 0) data can be fitted to a power-law, 
with TDOS (R = 0,w = 0) = a 0 \W - 1F/“| /3 . The obtained 
critical exponent for large enough clusters f3 ~ 1.62 ± 0.10 
is in good agreement with exact results. Note that the 
— f r(K, uj)dK.dui data for N c = 64 and 125 has been nor¬ 
malized with that of N c = 1. 


indicating no tendency towards localization with increas¬ 
ing disorder. Our results show that as N c increases, the 
critical disorder strength W c in the TMDCA system¬ 
atically increases until it converges to the exact value 
W c « 2.10iMi-i6 i 50-^2 i 60 at cluster size N c > 12. The 
cluster N c = 12 is the first cluster with a complete 
nearest-neighbor shell based on Betts cluster classifica¬ 
tion*^- From this cluster onward, W c converges to ~ 2.10, 
but the trajectory of the mobility edge in the re-entrance 
regime continues to change until it also converges to the 
exact results at larger N c . This effect is due to the sys¬ 
tematic incorporation of coherent backscattering as the 
cluster size increases and will be elaborated in more de¬ 
tail later. 

To extract the order parameter critical exponent (/3), 
we fit our data in Fig. [6] for the largest system size 
considered here ( N c = 125) using the power law: 
TDOS(w = 0) = ao\W — W/ lt |^. Following the proce¬ 
dure as explained in Sec. IIII El we obtain (3 ~ 1.62 ±0.10 
with a corresponding critical disorder strength from the 
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fit Wf lt ~ 2.23 ± 0.10. The fit overestimates the criti¬ 
cal disorder strength, as compared to the computed one, 
due to the difficulty in determining the scaling regime, 
as discussed below. This also causes the error bars to be 
larger than obtained from other methods. Nevertheless, 
the critical parameters from the fit are in good agreement 
with the recently reported value of /? ss 1.674& It is also 
in general agreement with the values listed in Table IIIII 

Apart from the typical density of states, the localiza¬ 
tion transition in the gapless single-particle excitations of 
the Anderson insulator can be studied using the return 
probability of an electron to a sitei 31 ' 58 The probability 
of quantum diffusion (or the return probability) describes 
the probability of a quantum particle (or a wave) to go 
from site l to l 1 in a time t. After disorder averaging, the 
return probability has basically three key contributions: 
a) the probability of going from site l to V without any 
scattering, b) the probability of going from site l to l' by 
an incoherent sequence of multiple scattering (known as 
diffusion) and c) the probability to go from site l to V via 
a coherent (or enhanced) multiple scattering processes 
(e.g., the cooperon ). 

In the DCA or TMDCA, it is more convenient to 
measure the probability of an electron remaining on a 



FIG. 7. (Color online). Right panels: The prob¬ 
ability that an electron remains localized at all times 
P(oo, 77 ) for the Anderson model for N c = 1 and 64 
at varying disorder strengths. We used the fact that 
Pit ^oo, V )= lim ^ E; fZ todu>'( A W,iW ), where 

A(l,uj) = — l/nQG(l,l,cu) is the local coarse-grained (but not 
disorder averaged) spectral function. As 77 0, P(oo,r/) ex¬ 

trapolates to zero for small disorder strength indicating metal¬ 
lic behavior but does not extrapolate to zero anymore as the 
disorder strength is systematically increased towards the crit¬ 
ical disorder strength leading to the transition (see the inset 
where this is manifestly illustrated). Left panels: The prob¬ 
ability of an electron on a site thereby remaining trapped at 
finite time t for N c = 1 and 64 for the same parameter as 
P(oo, 77 ) on a semi-log plot. 


given site l for all time P(t —> 00,77) = lim (\G(l, l, t)\ 2 ) 
= lim ^ fZ de(\G(l, l, e + ir/)l 2 )2& This will depend on 

77—^0 n J 00 ' ' 

the localization length, but if a significant fraction of 
the eigenstates of the disordered thermodynamic spec¬ 
trum are localized states, P(oo,? 7 ) is expected to be 
nonzero i 28 ! 31 ’ 53 : 59 Since the cluster is formed by coarse- 
graining the real lattice problem in K-space, there is a 
one-to-one correspondence between local quantities on 
the cluster and real spaced In Fig. [7] right panel, we 
show the P( 00 , 77 ) for the cluster sizes N c = 1 and 64 for 
various disorder strengths. As it is evident from the plot, 
for relatively small disorder strength W ~ 0.5, P(oo,r/) 
extrapolates to zero and becomes nonzero as the local¬ 
ization transition is approached. Just like the P(oo,?y), 
the finite time probability that an electron on a site l re¬ 
mains after some time t denoted as: P[t) = (\G(l, l, t)\ 2 ) 
is a vital parameter for detecting the localization of elec¬ 
trons. As shown in Fig. [7] for N c = 1 and 64 clusters, 
P(oo,rj) and P{t) contains the same information of the 
excitation spectra. In Fig. [3 left panel, we show the P(t) 
for the same parameters as P( 00 , 77 ). Hence, a character¬ 
istic finite long time P(t) denotes localized eigenstates. 
Again, systematic transition from a metallic regime (for 
small disorder) to an insulating regime (for a disorder 
strengths close to the critical value of the Anderson lo¬ 
calization transition in N c = 1 and 64 size clusters, re¬ 
spectively) is observed. Unlike in the DCA, the localiza¬ 
tion transition manifests clearly in the P(oo, 77) and P(t) 
calculated in the TMDCA since, even though the density 
of states (ADOS) calculated within the TMDCA is finite 
as in the DCA (the ADOS is a conserving quantity), the 
hybridization rate at the same energy depends highly on 
the strength of the disorder (cf. inset of Fig. [ 6 ]). In fact, 
it vanishes continuously with the disorder strength and 
goes to zero at the same point where the typical density 
of states vanishes. Hence, since the hybridization rate be¬ 
tween the cluster and the host vanishes continuously as 
the critical disorder strength is approached, the TMDCA 
method is able to capture the localization transition even 
when the ADOS calculated with the TMDCA is finite. 

The probability distribution function (PDF) is another 
natural quantity to characterize the 3D Anderson local¬ 
ization transition due to the fact that the “typical” value 
of a “random” variable corresponds to the most probable 
value of the PDFj 33 i 50 A proper description of electron 
localization in disordered systems requires consideration 
of the distribution functions for the quantities of inter¬ 
est^, so we calculate the PDF of the cluster-momentum- 
resolved DOS p(K,uj = Ik) (at different momenta cells 
K and energy u> = Ik) sampled over a large number of 
disorder configurations. 

In Figure [HI we show the evolution of the 
PDF[p(i7, ui = Ik)] with W. As is evident from the plot, 
for a relatively small disorder, the cells show a Gaus¬ 
sian distribution which gradually becomes log-normal 
and highly skewed as the critical disorder strength is ap¬ 
proached. 
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p(K,Q) = E ) 


p(K,a> = 


p(K,co = ej 


□—s N = 1 

c 



FIG. 8. (Color online). The evolution with disorder strength 
of the probability distribution of the density of states at differ¬ 
ent cluster cells for N c = 38. Utilizing the irreducible wedge 
property and particle-hole symmetry, the original 38 cells are 
reduced to 4 cells. For small disorder strength, the cells show 
Gaussian distributions whereas close to the critical disorder 
strength ~ 2.0, all the cells show log-normal distributions. 

The analysis of the spectra properties of the Ander¬ 
son mode l 39 ' 62 shows that for relatively small disorder 
strengths the states are still delocalized, and the am¬ 
plitude of the wave functions associated with them is 
more or less the same on every site. The distribution of 
the local DOS with respect to disorder configurations is 
Gaussian with the most probable value coinciding with 
the arithmetic mean average value. However, for suf¬ 
ficiently large disorder strength or in the proximity of 
the band tails, the spectrum consists mainly of discrete 
eigenvalues, and the associated eigenfunctions are expo¬ 
nentially localized with substantial weight only on a few 
sites. The distribution is therefore extremely asymmetric 
(log-normal), with a most probable value much smaller 
than the arithmetic mean value. At this point, most of 
the weight is concentrated around zero. As is evident 
from Fig. [51 we indeed observe such behavior in our re¬ 
sults. 

We show in Fig. [9] the phase diagram of the Ander¬ 
son localization transition in the disorder-frequency (W- 
w) plane constructed from our TMDCA procedure for 
the box disorder distribution. Here, we show the mobil¬ 
ity edge trajectories given by the frequencies where the 
TDOS vanishes at a given disorder strength W, and the 
band edge determined by the vanishing of the ADOS cal¬ 
culated within the DCA. 

As is evident from Fig. [HI at N c = 64 the W c at oj = 0 
is the same as that for N c = 125 but different from the 
N c = 1 case. This shows that the W c converges to ~ 2.10, 
while, the trajectory of the mobility edge continues to 
change with N c . This may be understood from the dif¬ 
ferent localization mechanisms for states at the band cen¬ 
ter and the edge i 51 ' 63 ' 64 Hence, for large enough clusters, 


FIG. 9. (Color online). The phase diagram of the Anderson 
localization transition in 3D for the box disorder distribution 
obtained from TMDCA simulations. A systematic improve¬ 
ment of the trajectories of the mobility edge is achieved as the 
cluster size increases. At large enough N c and within compu¬ 
tation error, our results converge to those determined by the 
transfer matrix method (TMM). The TMM data is for system 
sizes of length L = lx 10 5 (number of multiplications), the 
range of system widths used is M = [2,16] and reorthogo- 
nalization is done every 5 transfer matrix multiplications (see 
Appendix [A] for details). The error was determined in a point 
in the region of most discrepancy between the methods (up¬ 
per bound of re-entrance) by a finite size scaling analysis (see 
Appendix^. The black line with filled circles denotes the Lif- 
shitz boundaries (extracted from the ADOS calculated within 
the DCA). 


we are able to converge to the exact result. In partic¬ 
ular, as the cluster size N c increases, the mobility edge 
trajectories are systematically reproduced, with the re¬ 
entrance behavior gradually captured with large cluster 
sizes. As we increase the cluster size the DOS systemat¬ 
ically acquires states in the band tails, which are zero in 
the iV c =l case (it is well known that single-site theories 
like CPA or TMT do not capture such states). Accord¬ 
ing to Bulka et a/.; 47 ' 51 deep trapped states dictate the 
physics at large energies. Hence, by making N c > 1 we 
systematically inject additional states that tend to push 
the localization edge outward. States at the band center 
become localized mainly due to coherent backscattering 
while those above and below the bare band edges are 
initially localized in deeply trapped states. They become 
delocalized with increasing W due to the increasing DOS 
at these energies and hence increasing quantum tunneling 
between the deeply trapped states. They finally become 
localized again with increasing disorder, which explains 
the re-entrant behavior. Since coherent backscattering 
requires a retracing of the electronic path, the effective 
length scales captured by the cluster are doubled, such 
that W c converges very quickly at the band center. On 
the other hand, the quantum tunneling mechanism has 
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FIG. 10. (Color online). The scaling of the imaginary typical 
hybridization function (— £yr(K,w)) for the 64 site cluster at 
disorder strengths of W = 1.8 and 2.0 on a semi-log scale. 
The labels A-F and their associated momenta K correspond 
to each of the six distinct cells obtained using the point-group 
and particle-hole symmetry (p( K,w) = p(Q — K, — u), with 
Q = (tt, 7r, 7 t)) of the cluster). Observe that the mobility 
edges may be collapsed on top of each other by multiplying 
each of the hybridization functions with a constant such that 
—AF(K,aj) = —a(K) x $jT(K,u;), where a(K) is a scaling 
constant, in agreement with Mott’s idea of energy selective 
Anderson localization transition. 


no path doubling and requires multiple deeply trapped 
states in the cluster and, therefore, converges more slowly 
with N c . 

Figure [TO] is the scaling of the imaginary part of the 
typical hybridization function (—QF(K,u;)) for N c = 64 
at W = 1.8 and 2.0, scaled by the factor a(K), so that the 
tails overlap. Even though different K-cells go to zero at 
different rates, they each share the same unique mobility 
edge. Since the local disorder potential induces elastic 
scattering of a state at any momenta into any other with 
the same energy, there will be mixing of the localized 
and the extended states at the same energy. As a result, 
the localized and extended states cannot coexist at the 
same energy. Then the mobility edge can only exist at 
the point where all the states in each cell on the cluster 
are localized. As shown in Fig. [10] the collapse of the 
tails for all K such that -QT(K,w) = -a(K)3T(K, w), 
where a(K) is a scaling constant, validates Mott’s idea 
of energy selective Anderson localization! 65 ' 66 


B. Alloy Model 

The application of the disordered tight-binding Hamil¬ 
tonian © to the alloy model represents one of the most 
studied physical systems. This stems from the fact that 
the two potentials energies W a and W& depict the poten¬ 
tial landscape of e.g., a binary alloy A Ca B\_ Cal with each 
of the sites occupied either by atom “A” or “B” with 


concentrations c a and q, = 1 — c a , respectively. 



FIG. 11. (Color online). Comparison of the average (typical) 
density of states calculated with the DCA (TMDCA) and 
the kernel polynomial method (KPM) for the binary alloy 
system (144 = —144) at various values of the local potential 
144 and concentrations c a for the cluster size A4=125. The 
kernel polynomial method uses 2048 moments on a 48 3 cubic 
lattice, and 200 independent realizations generated with 32 
sites randomly sampled for each realization. 

To explore the applicability of our method for the 
study of binary alloys, we start the discussion of this 
section by showing in Fig. |TT] the calculated typical (av¬ 
erage) density of states from the TMDCA (DCA) pro¬ 
cedure as compared to the TDOS (ADOS) calculated 
within the for various concentrations and 

disorder strengths. The importance of the TDOS is ev¬ 
ident since for all the disorder strengths and concentra¬ 
tions, the ADOS remains finite around the two energies 
W a and Wb, while the TDOS at a fixed concentration 
vanishes continuously with the strength of the disorder 
with smaller values in the sub-band with the lowest con¬ 
centration. Fixing the strength of the disorder and vary¬ 
ing the concentration, the sub-bands with the smallest 
concentration have fewer states. We note that there are 
subtle differences between the results for N c = 64 (finite 
cluster) and single-site N c = 1 (CPA) (not shown) due to 
the incorporation of spatial correlations in the finite clus¬ 
ter which are missing in the local CPA. In fact, the TMT 
underestimates the extended region and misses small but 
important nonlocal features in the spectra.— 

To further benchmark our results for the binary al¬ 
loy model, let us focus on the comparison of the average 
(typical) DOS calculated with the DCA and TMDCA 
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• TMM 
— TMT (N c =l) 



FIG. 12. (Color online). Disorder-energy (W a -w) phase di¬ 
agram of the Anderson localization transition in 3D for the 
binary alloy system A Ca Bi- Ca for c a = 0.5. Observe the sys¬ 
tematic improvement of the trajectories of the phase diagram 
for the clusters (N c = 64 and 216) in basic agreement with 
the numerical results from TMM. The system widths used for 
TMM are M=[6,12] and the length is scaled with the width 
as L = M x 10 4 . See Appendix lAl for details. 

(N c = 64) with the KPM data. As can be seen in Fig, fill 
the TMDCA and DCA results reproduced those from the 
KPM, showing that our formalism offers a systematic 
way of studying the Anderson localization transition in 
binary alloy systems. Such a remarkable agreement is an 
indication of a successful benchmarking of the TMDCA 
method^ 

We extract the mobility edges shown in Fig. [l^] by find¬ 
ing the energy where the TDOS vanishes at a given value 
of the disorder potential. As can be seen in Fig. [12] the lo¬ 
cal N c = 1 boundaries are narrower than those obtained 
for the finite cluster indicating that the TMT strongly 
underestimates the extended state regime. The compar¬ 
ison of the mobility edge boundaries for the TMDCA 
with those obtained from the TMM calculations show 
quite good agreement. This again is a confirmation of a 
successful benchmarking of the TMDCA for treating the 
binary alloy model. At the center of the band and for 
c a = 0.5, we obtain a critical disorder strength W c of « 
0.7 for the TMDCA in good agreement with the TMM 
(cf. Fig. fl2l). 

C. Gaussian Disorder Distribution 

The Gaussian (or normal) (Eq. I2bl) is a unique dis¬ 
tribution that other disorder distributions are built on 
and has many physical applications including the study 
of molecular-doped polymersi 68 ’ 69 

To further explore the versatility of our method, we 
apply it to study systems with the disorder defined by 



FIG. 13. (Color online). Comparison of the average (typ¬ 
ical) density of states calculated with the DCA (TMDCA) 
and the kernel polynomial method (KPM) for the Gaussian 
disorder distribution at various values of the disorder strength 
for cluster size N c = 64. The kernel polynomial method uses 
4096 moments on a 48 3 cubic lattice, and 1000 independent 
realizations generated with 32 sites randomly sampled from 
each realization. 

the Gaussian distribution function (Eq. I2bl) . Again, we 
use the typical density of states (TDOS) as the order pa¬ 
rameter and the transition to the Anderson insulator is 
obtained at the disorder strength where the TDOS van¬ 
ishes. The typical and average density of states obtained 
from the TMDCA and DCA, respectively, and those ob¬ 
tained from the kernel polynomial method are shown in 
Fig. [13] for various values of the disorder strength. As 
can be seen, the TDOS at all frequencies systematically 
goes to zero as the disorder strength increases while the 
ADOS remains finite. Again, our TMDCA formalism re¬ 
produces accurately the results from the kernel polyno¬ 
mial method. We note some subtle differences between 
the TDOS calculated from the TMDCA and the KPM 
while there are no noticeable differences in the average 
density of states from the DCA and KPM. This may 
be due to the finite broadening utilized in the KPM, 
which contributes additional tails to the already expo¬ 
nential tails of the TDOS. We remark that aside from 
the small initial broadening value (~ —0.01) used in the 
initialization of the TMDCA at the very first iteration, 
no broadening parameter is utilized for later iterations. 

We show in Fig. [TT] the evolution of the typical den¬ 
sity of states TDOS(w = 0) at the band center as a func¬ 
tion of disorder strength for the local TMT (N c = 1) 
and the TMDCA ( N c = 64 and 125). Our results in¬ 
dicate that the critical disorder strength (defined as the 
W where the TDOS vanishes) systematically increases 
as the cluster size is increased converging to W c ~ 5.30 
as soon as the size of the cluster N c > 12. This is 
in good agreement with the numerically exact values of 
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FIG. 14. (Color online). The TDOS(u; = 0) vs. disorder 
strength W at different cluster sizes N c = 1,64, and 125 for 
the Gaussian disorder distribution. The TDOS(u> = 0) van¬ 
ishes at the critical disorder strength W c - At N c = 1, the 
critical disorder strength Wv° =1 ~ 4.0. As the cluster size 
increases, the critical strength systematically increases with 
yj/jv c »i 2 ^ g 3 q o.lO, showing a quick convergence with 
the cluster size. The inset shows the integrated hybridiza¬ 
tion function (—3 f T(K, uj)dKduj) as a function of disorder 
strength W. Observe that —51 / T(K, Lj)dK.du> vanishes at the 
same disorder strength as the TDOS. The dashed lines are the 
—SS f T(K, uj)dK.duj from the DCA. Observe that it is a con¬ 
stant regardless of the disorder strength. This shows that the 
DCA, even though it incorporates spatial correlations does 
not describe the Anderson localization transition. Moreover, 
near the critical region the TDOS(cj = 0) data can be fitted to 
the power-law TDOS(w = 0) = ao |W — I' 3 . The obtained 

critical exponent for large enough clusters /3 « 1.57 ± 0.10 is 
in good agreement with the numerically exact results.— Note, 
the —Sr f T(K, uj)dKdu> data for N c = 64 and 125 has been 
normalized with that of N c = 1. 



FIG. 15. (Color online). The phase diagram of the Ander¬ 
son transition in 3D for the Gaussian disorder distribution 
obtained from TMDCA simulations. A systematic improve¬ 
ment of the trajectories of the mobility edge is achieved as 
the cluster size increases. At large enough N c , our results 
converge to those determined by the TMM which are calcu¬ 
lated for widths M = [2,16] and for a system length (number 
of transfer matrix multiplications) L = lx 10 5 , where the 
matrix products are reorthogonalized every 5 transfer matrix 
multiplications (see Appendix El. The deviation between the 
TMDCA and TMM results is consistent with the behavior 
observed for the box and Lorentzian disorders and can be at¬ 
tributed to the fact that a finite grid in energy is used for the 
TMDCA which tends to cause the typical density of states 
to be larger, hence slightly overemphasizing the metallic be¬ 
havior and as such, the mobility edge is slightly larger when 
compared to TMM in certain frequency ranges near the band 
edge. The effect is most pronounced here due to the small 
density of states near the band edge. In addition, near the 
re-entrance regime, the TMM also has difficulties due to an 
increase in finite size effects (see Appendix 


5.225±0.125^£~— and 5.32— Fitting our data for the 
largest system size considered here ( N c = 125) using the 
power law: TDOS(w = 0) = clq\W — Wl lt \P and following 
the procedure as explained in Subsection IIII El we obtain 
the order parameter critical exponent /3 ~ 1.57 ± 0.10 
with a corresponding critical disorder strength from the 
fit of W/ lt ~ 5.53 ± 0.10. This value of /? is in good 
agreement with the value we obtained for the uniform 
disorder distribution (cf. Tabic [TT| and in good agree¬ 
ment with the recently reported value of /3 ~ 1.67— It 
is also in general agreement with the values listed in Ta¬ 
ble Ell The good agreement between the j3 we obtained 
from the uniform and Gaussian disorder distribution is 
a manifestation of the universal nature of the Anderson 
localization transition— iZS 

To explore the trajectories of the mobility edge for the 
Gaussian disorder distribution, we show in Fig. [15] the 
phase diagram in the energy-disorder plane for various 
cluster sizes as compared to the TMM result— For any 
given disorder strength W, the mobility edge is defined by 
the frequency where the TDOS vanishes. Unlike the crit¬ 


ical disorder strength which converges quickly with the 
cluster size N c > 12, the trajectory of the mobility edge 
continues to change with N c converging to almost the nu¬ 
merically exact results for N c = 125. The physical rea¬ 
sons for the quick convergence of W c and the progressive 
change of the mobility edge with the cluster size are the 
same as those described above for the box disorder distri¬ 
bution (cf. Section IIII Al ). As can be seen, the single-site 
TMT underestimates the extended region just as in the 
previously presented disorder distributions. There are, 
however, some subtle differences between our data and 
the TMM results around the re-entrance regime. The 
cause of this difference will be discussed in Section IIII FI 

D. Lorentzian Disorder Distribution 

We next apply our TMDCA formalism to study sys¬ 
tems with the Lorentzian (or Cauchy) disorder distribu¬ 
tion (cf. Eq. I2c|) . We show in Fig. [Tn] how the band 
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FIG. 16. (Color online). The TDOS(u; = 0) vs. disorder 
strength W at cluster sizes JV C = 1, 64, and 125 for the 
Lorentzian disorder distribution. The TDOS(a> = 0) van¬ 
ishes at the critical disorder strength W c - At N c = 1, the 
critical disorder strength W^ c=1 fts 0.6. As the cluster size 
increases, the critical strength systematically increases until 
He >12 ~ 0.94 ± 0.10, showing a quick convergence with 
the cluster size. The inset shows the integrated hybridiza¬ 
tion function (—f T(K, uj)dKduj) as a function of disorder 
strength W. Observe that —£r f T(K, uj)dK.du> vanishes at 
the same disorder strength as the TDOS. The dashed lines 
are the —9 f r(K,cj)dKdw from the DCA. Observe that it 
is a constant regardless of the disorder strength in agreement 
with the observations we have for other disorder distributions. 
We can fit the TDOS(cu = 0) data near the critical region 
to the power-law, with TDOS(w = 0) = ao\W — 

The obtained critical exponent for large enough clusters 
P « 1.60 ± 0.10 is in good agreement with numerically ex¬ 
act results4£ Note, the —O' f T(K, uj)dK.du> data for N c = 64 
and 125 have been normalized with that of N c = 1. 


center of the typical density of states (TDOS(u> = 0)) 
changes as the disorder strength is increased for the lo¬ 
cal TMT (N c = 1) and the TMDCA (N c = 64 and 
125). As can be seen from Fig. [16j our results de¬ 
pict that the critical disorder strength systematically 
increases as the cluster size increases converging to 
W c ~ 0.94 for N c > 12 in good agreement with the 
numerically exact values of 0.95±0.125^— and 1.07^ 
respectively. Fitting our data for the largest system 
size considered here (N c = 125) using the power law: 
TDOS(w = 0) = ao\W — W ^ lt (see Subsection IIII El 
for the description of how the (3 is extracted), we can in¬ 
fer the order parameter critical exponent /? ~ 1.60±0.10 
with a corresponding critical disorder strength from the 
fit, Wl lt ~ 0.97 ± 0.10. The obtained /3 from our fit 
is in good agreement with the values listed in Table Hill 
and in good agreement with the recently reported value 
of ft « 1.674^ This /3 value of the Lorentzian disorder 
distribution is also in good agreement with the values we 
obtained for the box and Gaussian disorder distributions, 
respectively. This further illustrates the universal nature 
of the Anderson localization transition ; 49 ! 70 



FIG. 17. (Color online). The phase diagram of the Ander¬ 
son localization transition in 3D for the Lorentzian disorder 
distribution obtained from TMDCA simulations. The TMM 
data were computed for a system length of L = 8 x 10 5 for 
system widths of M = [2,16] and re-orthogonalization is done 
every four multiplications (see Appendix [A] for details). The 
discrepancy between the TMDCA and TMM is due to the na¬ 
ture of the Lorentzian distribution as explained in Sec. IIII FI 


We conclude our study of the application of the typ¬ 
ical medium dynamical cluster approximation to the 
Lorentzian disorder distribution by presenting in Fig. [T7] 
the phase diagram in the energy-disorder plane. Unlike 
the box and Gaussian disorder distributions, our simu¬ 
lations show that the Lorentzian distribution does not 
have re-entrance of the mobility edge. The lack of re¬ 
entrance of the mobility edge in the Lorentzian disorder 
distribution may be attributed to the absence of finite 
variance in this form of distribution. For the single-site 
CPA (N c = 1), the critical parameters are woefully un¬ 
derestimated. However, we systematically converge to 
the numerically exact results as the size of the cluster is 
increased. As it is obvious from Fig. [T71 for as small as 
N c = 64, we converge almost to the exact TMM results^ 
We again see remarkably good agreement between our ef¬ 
fective mean-field method for the Anderson localization 
transition and the numerically exact calculations. 


E. Critical Parameters 

The critical parameters, including the critical disorder 
strength W c and the order parameter critical exponent 
P are summarized in Table EH for different cluster sizes 
using the box disorder distribution as a case study. W^ al 
was determined as the W where the TDOS (w = 0) van¬ 
ishes. Observe that as N c increases, lb” 1 systematically 
increases with W” 1 « 2.10 ± 0.01, showing a quick con¬ 
vergence with N c . 

The order parameter critical exponent (j is obtained by 
fitting the power law: TDOS(w = 0) = ao x \W — 
directly to our data and systematically searching for 
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TABLE II. The calculated and fitted critical disorder 
strengths W^ al and W/ u and the order parameter critical 
exponent /? obtained from our fit for the box disorder distri¬ 
bution. W^ al is defined by the vanishing of the TDOS (ui = 0). 
/3 and W/ 1 * are obtained by fitting the TDOS(a> = 0) data 
with a power law, TDOS(a) = 0) = ao|VF — W/ it \^. 


N c 

w t al 

W f c u 

P 

1 

1.66T0.01 

1.65T0.10 

0.96T0.10 

64 

2.10T0.01 

2.18T0.10 

1.46T0.10 

125 

2.10T0.01 

2.23T0.10 

1.62T0.10 


the best data point away from the transition where 
the fit still follows the actual data (scaling regime). 
This becomes imperative since away from the tran¬ 
sition, the data are not expected to fit the form 
TDOS(w = 0) = oo x |W — W'f lt \^ and also close to the 
transition, there should be a crossover to a mean-field 
form since the TMDCA treats the longest length scales 
in a mean-field approximation. So, the fit may only be 
done between these limits. The ambiguity in the deter¬ 
mination of the fitting region increases our error bars on 
the exponent, and causes W/ lt to be overestimated. In 
addition, the strong fluctuation of the TDOS in the prox¬ 
imity of the critical point also increases our error bars. 
We show in Fig. [18] a comparison plot of the fit and our 
data for the 125 site cluster for the box, Lorentzian, and 
Gaussian disorder distributions, respectively. The fit of 
the power law to the scaling region of the data gives the 
value of the ft in an unambiguous manner. The obtained 
values of /? from the various cluster sizes, for instance, for 
the box disorder distribution are shown in Table [TT] One 
can see that our /3 systematically approaches the numer¬ 
ical experimental valu o 44 ’ 46 ^ 50 for large enough clusters 
(here, largest N c simulated is 125) as listed in Table Hill 
for the box, Gaussian, and Lorentzian disorder distribu¬ 
tions, respectively, and in comparison with other numer¬ 
ical values. 


F. Difficulties in Extracting the Mobility Edge at 
Higher Disorder 

As explained in the previous sections, the mobility edge 
is obtained by locating the frequencies where the typical 
density of states vanishes at a given disorder strength 
W. The band edge is determined by the vanishing of 
the average density of states calculated within the DC A. 
As it is obvious from the phase diagrams for the various 
disorder distributions CFigs. Ii)l [T5l andlTTl respectively), 
there are some discrepancies between the phase diagram 
obtained within the typical medium dynamical cluster 
approximation and the transfer matrix method. They 
can be attributed in part to the form of the disorder dis¬ 
tribution. For example, the bare DOS for the Gaussian 
and Lorentzian disorder distributions are known to have 
exponential tails. The severity of the exponential nature 


of the tails associated with the various disorder distribu¬ 
tions increases as box, Gaussian, and Lorentzian, in that 
order. At these higher disorder strengths (once a mobil¬ 
ity edge develops), the TDOS naturally develops tails. 
We note that in our computations, apart from the initial 
small broadening value ~ —0.01 used in the initialization 
of the self-energy (needed only for the first iteration), no 



FIG. 18. (Color online). The band center of the typical 
density of states (TDOS(f? = 0; to = 0)) at various dis¬ 
order strengths (W) for the 125 site cluster for the Box, 
Gaussian, and Lorentzian disorder distributions. The lin¬ 
ear region of the data can be fit with a scaling ansatz: 
TDOS(w = 0) = a 0 x \W - with /3 ~ 1.62 for the 

Box disorder distribution, p ~ 1.57 for the Gaussian disorder 
distribution, and /5 ~ 1.60 for the Lorentzian disorder distri¬ 
bution in good agreement with the recently reported value of 
/? ~ 1.67.—^ The circles in the plot depict the data point 
where the fit starts to deviate from the data. 



FIG. 19. (Color online). The plot of the typical density of 
states (TDOS) at uj = 0.0 and w = 1.3 for various disorder 
strengths for the (a) Box, (b) Gaussian, and (c) Lorentzian 
disorder distributions. 















17 


TABLE III. Critical parameters of the Anderson localization for the various studied disorder distributions in 3D obtained using 
the TMDCA in comparison with the numerically exact results. We use 4t = 1 as our energy unit. We note that the critical 
exponents (/3 and u) are independent of disorder distribution (universal) as verified by the multifractal analysis,— the analytic 
results of Wegner^S and detailed finite size scaling— Abbreviations used in the Table are: transfer matrix method (TMM), 
multifractal finite size scaling (MFSS), level statistics (LS), kicked rotor (KR), and experimental atomic waves (Exp-AW). 


Author 

Critical Disorder 

Critical Exponent 

Method 


Bo Ga 

Lo 

V 

0 


Present study 

2.10T0.10 5.30T0.10 

0.94±0.10 

- 

1.57-1.62 

TMDCA 

Slevin et al 

2.067 

1.067 

1.573-1.577 

- 

TMM 

Slevin et aF 

2.068-2.073 5.32 

1.066 

1.58 

- 

TMM 

Slevin et al / 

2.056T0.014 

- 

1.59-1.60 

- 

TMM 

Rodriguez ^ 

2.066-2.067 

- 

1.59 

1.65-1.68 

MFSS 

Rodriguez ^ 

2.066-2.071 

- 

1:58±0.03 

- 

MFSS 

Markos^ 

2.063,2.067 

- 

1.47-1.55 

- 

MFSS 

MacKinnon^ 

2.063 ±0.05 

- 

1.54±0.08 

- 

TMM 

MacKinnon et al? 

2.063 ±0.063 

- 

1.2±0.3 

- 

TMM 

Bulka et al._ 

2.038-2.063 5.23±0.13 

0.95±0.13 

- 

- 

TMM 

Milde et al™. 

- 

- 

1.62±0:07 

- 

TMM 

Shklovskii et al A 

2.0±0.063 

- 

1.50±0.15 


LS 

Zharekeshev et aZT 

2.05 

- 

1.4±0.15 

- 

LS 

Hofstetter ^ 

2.719±0.012 

- 

1.35±0.15 

- 

LS 

Lopez ^ 

- 

- 

1.63±0.05 

- 

KR 

Grussbach et al 7 

2.02 5.23 

- 

1.32±0.02-1.37±0.02 

1.32±0.02-1.37±0.02 

MFSS ^ 

Lemarie ^ 

- 

- 

1.58±0.01-1.60±0.03 

- 

KR^ 

Lemarie ^ 

- 

- 

1.40±0.30 

- 

Exp-AW ^ 


The order parameter critical exponent /3 and the correlation length critical exponent v can be transformed from one to the 
other using the hyperscaling relation of Ref. (0) /3 = (q„ — d)v, where a 0 is the Lipschitz-Holder exponent, which gives the 
maximum value of the multifractal spectrum. The most recent estimates as reported in Ref. (0) are a 0 = 4.048, v = 1.59 
and /3 = 1.67. 


a The authors obtain a 0 = 4.0 such that in the hyperscaling relation /3 = (a 0 — d)v, 0 equals v. 

b Authors of this paper show that their quasiperiodic kicked rotor belongs to the same (orthogonal) universality class as the ‘random’ 
Anderson model. 

c Authors of this paper reported that experiments were done on the atomic kicked rotor by a sequence of kicks to the atomic cloud and 
measure its dynamics. 



broadening factor is utilized. As such, these tails that 
emerge as the disorder strength is increased towards W c 
are physical tails since the top and bottom of the bands 
will localize first. To demonstrate this, we show in Fig. [T9] 
a plot of the TDOS at uj = 0.0 and w = 1.3 for the box, 
Gaussian, and Lorentzian disorder distributions, respec¬ 
tively. Note the u = 1.3 frequency is arbitrary but chosen 


such that it is close to the re-entrance region of the mo¬ 
bility edge. As can be seen from the plots, for the box 
disorder distribution, even though there are small tails, 
the TDOS at w = 0.0 and w = 1.3 behave alike and dif¬ 
fer only in magnitude. There are no obvious long tails 
in either frequency ( u> = 0.0 and 1.3) that may mask the 
detection of the position of the mobility edge energies. 
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However, for the Gaussian and Lorentzian disorder dis¬ 
tributions, there are exponentially long tails especially at 
w = 1.3 that make pinpointing the exact position of the 
mobility edge energies highly non-trivial. This should 
further be understood from the fact that unlike the box 
disorder distribution, the Lorentzian disorder distribu¬ 
tion naturally has tails that decay very slowly at infinity, 
as |a:| , such that aside the zeroth moment (the area 

under the curve), all other higher moments do not even 
exists. 

We note that this difficulty is generic not only to the 
TMDCA but also for any method where the extraction 
of the trajectories of the mobility edge is based on the 
TDOS. A notable example is in the kernel polynomial 
method, which even shows more severe discrepancies (not 
shown). 

This deviation between the TMDCA and TMM at 
higher disorder strengths can also be attributed to the 
fact that the TMDCA, just like the KPM utilizes a finite 
frequency grid (which biases more towards the metallic 
regime) in contrast to the TMM which calculates the 
transmission of electrons at fixed energies. Even so, we 
note that the TMM also has its own shortcomings away 
from the band center due to the strong fluctuations in 
the Kramer-MacKinnon scaling parameter A (see Ap¬ 
pendix [A] for its definition) as can be seen from the phase 
diagrams for the various disorder strengths around the 
re-entrance regime. 


IV. CONCLUSIONS 

Here, we present a detailed study of the Anderson lo¬ 
calization transition using the recently developed typical 
medium dynamical cluster approximation (TMDCA) for 
the box, Gaussian, Lorentzian, and binary alloy disor¬ 
der distributions in three dimensions. For each distri¬ 
bution we find the TMDCA to be a successful, causal, 
numerically efficient, self-consistent and rapidly conver¬ 
gent method for the study of localization in disordered 
electron systems. 

With our formalism, we demonstrate that the typical 
DOS vanishes for localized states and is finite for ex¬ 
tended ones. Employing the typical DOS as an order pa¬ 
rameter for Anderson localization, we have constructed 
the disorder-energy phase diagram, extracted the order 
parameter critical exponent (/3) for each disorder distri¬ 
bution, and benchmarked them in good agreement with 
other numerically exact methods. Within our precision, 
we find that /3 for the Anderson localization transition 
is a universal parameter independent of disorder distri¬ 
bution in agreement with the multifractal analysis.— For 
distributions with a finite variance (box and Gaussian), 
we demonstrate that there are extended states outside 
the unperturbed band. 

We further show using the DCA (which includes spatial 
correlations) and a variant of the typical medium theory 
(which includes spatial correlations but suffers from self¬ 


averaging), the importance of the effective medium to 
properly characterize the Anderson localization transi¬ 
tion. We also demonstrate the inability of the single-site 
CPA and the TMT methods to accurately capture the 
localization and disorder effects in both the average and 
the typical DOS. We note that the single-site TMT, while 
being able to qualitatively capture the localization tran¬ 
sition, strongly underestimates the extended regions and 
fails to capture the critical parameters including the mo¬ 
bility edge trajectories and the exponents. In contrast, 
the TMDCA captures nicely the trajectories of the mo¬ 
bility edge with great improvement in the critical order 
parameter exponent. Most significantly, the TMDCA re¬ 
sults are in a quantitative agreement with exact numeri¬ 
cal results. 

The TMDCA formalism is computationally inexpen¬ 
sive and straightforward to implement since it requires 
only the computer time needed to invert small clusters 
(e.g., N c = 1-125), average over the disorder configura¬ 
tions, and iterate to convergence. Since only a small clus¬ 
ter is needed to get reliable data, material specific details 
may be incorporated. Once combined with electronic 
structure calculations^ and more sophisticated many- 
body techniques for electron interactions, it will open a 
new avenue for studying the localization phenomenon in 
real materials as well as the competition between disorder 
and electron correlations. To demonstrate the high effi¬ 
ciency of the TMDCA, as only small clusters are needed 
to get a converged result in good agreement with the 
TMM data, we compare the relative CPU time needed for 
the largest system size simulated in the TMDCA and the 
TMM. For the largest cluster size used in the TMDCA 
calculations, which is N c = 216, the computation time is 
~4 hours (running on a single processor). While for the 
TMM, which is perfectly parallel in both disorder and 
frequency, each point in the phase diagram can require 
significant computational effort. For example, the sys¬ 
tem sizes in Fig. [9] used ~20 hours on 64 processors per 
frequency. Since a separate TMM calculation is needed 
for each frequency, achieving the energy resolution of a 
typical TMDCA calculation (for a certain number of fre¬ 
quency grid points) would require the product of the 
number of processors used to parallelize over disorder 
times the number of grid points. As a calculation of a 
real material would require even larger system sizes than 
used here in the TMM, the TMDCA would prove much 
more computationally efficient for the purpose of study¬ 
ing real materials. 
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Appendix A: Details of the Transfer Matrix Method 

We benchmark our TMDCA results for the mobil¬ 
ity edge with the transfer matrix method (TMM) 
which is an established numerical method for determining 
the mobility edge by computing the localization length 
in disordered quantum systems. TMM is based on an 
iterative formulation of the Schrodinger equation where 
the wave function amplitude is computed at each site 
in a quasi-one-dimensional “bar” of length L and width 
M by successive multiplications of the transfer matrix 
that describes the transmission between each “slice” of 
size M d_1 . Thus, the Lyapunov exponent that measures 
the exponential decay of the wavefunction is explicitly 
computed, yielding the localization length. The system 
length L is the total number of transfer matrix multipli¬ 
cations. The numerical instability of the repeated multi¬ 
plications is avoided by periodically orthogonalizing the 
transfer matrix product with a Lapack QR decomposi¬ 
tion after a finite number of iterations.— The transfer 
matrix method finite size effects are larger for weak dis¬ 
order where the states decay slowly with distance and so 
have large values of the localization length, which results 
in more pronounced fluctuations in the data. Notice that 
the CPA and the DCA (same as the TMDCA) do not suf¬ 
fer such finite size effects for small disorder and are, in 
fact, exact in this limit. 

The mobility edge is obtained by calculating the di¬ 
mensionless Kramer-MacKinnon scaling parameter Am, 


which is the localization length divided by M *2 \ M 
scales as a constant for M —)■ oo at the transition— Pre¬ 
cise values of the critical disorder may be measured di¬ 
rectly from the crossing plots of A m vs. W. A finite size 
scaling is performed by expanding Am near the critical 
point using 

A m = /(M/e) » A c +a l \W-W c \M 1 / v +a 2 \W-W c \ 2 M 2 /' J +- ■ ■ 

(Al) 

and the data is fit using a least squares procedure— 

The data used in Fig. [5]was a third order polynomial in 
|W — W C \M 1 ' V . The critical disorder strength W c from 
the finite size scaling is averaged over many generated 
data sets via a bootstrap procedure— Any errors quoted 
in the TMM data are from the difference in the measured 
critical disorder from the finite size scaling analysis^ and 
the critical disorder measured directly from the crossing 
plots of Am vs. W. 


Appendix B: Details of the Kernel Polynomial 
Method 

To further benchmark our results, we utilize the kernel 
polynomial method to calculate the local DOS—2 In the 
kernel polynomial method analysis, instead of diagonaliz¬ 
ing the Hamiltonian directly, the local DOS is expanded 
in terms of an infinite series of Chebyshev polynomials. 
In practice, the truncated series leads to Gibbs oscilla¬ 
tions. The kernel polynomial method damps these os¬ 
cillations by a modification of the expansion coefficients. 
We use the Jackson kernel following previous studies on 
the Anderson models 
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